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Abstract 

We study the Taylor expansion for the solution of a differential equation driven by a multi- 
dimensional Holder path with exponent (3 > 1/2. We derive a convergence criterion that enables 
us to write the solution as an infinite sum of iterated integrals on a nonempty interval. We 
apply our deterministic results to stochastic differential equations driven by fractional Brownian 
motions with Hurst parameter H > 1/2. We also prove that by using L 2 estimates of iterated 
integrals, the criterion and the speed of convergence for the stochastic Taylor expansion can be 
improved using Borel-Cantelli type arguments when H £ (i, ~j. 
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1 Introduction 



This paper is divided into two parts. In the first part, we consider a deterministic differential 
equation of the following type: 

( d 

dX t = Y,Vi(X t )dyi 

* i=0 \ • ) 

< X = x 

where the ViS are Lipschitz C°° vector fields on R n with Lipschitz derivatives and where the driving 
signal y : M> -> M d+1 is /3-H61der continuous with /3 > 1/2. Since y : M> -> M. d+1 is /3-H61der 

with /3 > 1/2, the integrals / Vi{X s )dy\ are understood in Young's sense (see [T2] and [T3]) and it 

J o 

is known that the above equation has a unique solution that can be constructed as a fixed point of 
the Picard iteration in a convenient Banach space of functions (see [S], |13| and |1U|). 
We express, under suitable conditions, the solution of the differential equation as an infinite series 
involving the iterated integrals of the Holder path. To study the question of the convergence of 
the infinite series, which we will call the Taylor expansion associated with the differential equation 
later on, we will develop some convergence criteria based on an explicit upper bound of the iterated 
integrals of the Holder path. The main tool we use to obtain such bounds is the fractional calculus 
and the main assumption we need to perform the expansion is the analycity of the vector fields 
ViS in a neighborhood of the initial condition xq. To illustrate our results, we consider examples of 
differential equations on a Lie group. 

The second part of this paper is concerned with a stochastic differential equation (in short SDE) of 
the following type: 

( d 

( dX t = Vi (X t ) dB\ 

8=0 
^X = Xq 

where the Vj's are Lipschitz C°° vector fields on M. n with Lipschitz derivatives and where (Bt)t>o 
is a d-dimensional fractional Brownian motion with Hurst parameter H > 1/2. Throughout the 
paper, we use the convention that B® = t. We will be interested in the convergence of Taylor 
expansion of the solution to the SDE. The deterministic estimates we obtained in the first part 
will automatically lead to convergence criteria for the stochastic Taylor expansion. The rate of 
convergence to and L p estimates of the remainder term are then studied. Afterwards, inspired 
by some ideas of Ben Arous |TJ (see also Castell [5j), we will also use probabilistic methods to 
obtain a convergence criterion, which improves the convergence radius of the stochastic Taylor 
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expansion in some cases. The key point to our probabilistic methods are the L2 estimates of the 
iterated integrals of the fractional Brownian motion that were obtained by Baudoin-Coutin in [3] 
and Neuenkirch-Nourdin-R6Bler-Tindel in [7]. 

2 Taylor expansion for differential equations driven by Holder 
paths 

In this section, we present our main theorem. We derive a convergence criterion, under which, we 
can, on a nonempty time interval, write the solution of our differential equation (jl.ip as an infinite 
series of iterated integrals of the Holder path y. The materials in this section will be organized in 
the following way: 

1. In Section 2.1, we will recall some basic facts about fractional calculus that will be needed in 
our analysis. 

2. In Section 2.2, we will briefly discuss the existence and uniqueness of the solution to the 
differential equation The main result is borrowed from [8]. We will also define, in 
apparently two different ways, the Taylor expansion of the solution of the differential equation 
(jl.ip . The two apparently different expansions are then shown to eventually be the same. 

3. In Section 2.3, we will state our main theorem. 

4. In Section 2.4 we will apply our results to the case of differential equations on Lie groups. 
2.1 Preliminaries : Fractional calculus 

Let us first do some remainders about fractional calculus. For further details, we refer the reader 
to ([I!]) or ([!]). Let / G L l (a,b) and a > 0. The left-sided and right-sided fractional integrals of 
/ of order a are defined by: 



respectively, where (—1) a = e t7va and T(a) = J °° u a 1 e u du is the Gamma function. Let us 
denote by (respectively I^_{LP)) the image of L p (a,b) by the operator I" + (respectively 

J2+). If / G Ia+( LP ) (respectively / G I£_(L P )) and < a < 1, we define for x G (a, b) the left and 
right Weyl derivatives by: 




and 
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1 / ff x ~\ r x f( x \ - fty) \ 

D ^ nx) = W^)V^^ + a !a {x-y)^ dv ) 1(a ' b){x) 

and respectively, 

n° ft \ (~l) a ( m , f b f(x)-f(y) , \ 

Now recall that from ([8]), for a parameter < a < 1/2, W T -a ' oo (0, T) is defined as the space of 
measurable function g : [0, T] — >■ M such that: 



Clearly, for every e > 0, 



\ g(t)-g(s)\ [* \g(v)-g (s)\ 



a 1 - a+e (o,T) c w^~ a '°°(o,T) c c 1 - a (o,r) 

Moreover, if 5 e H^~ Q,oo (0, T), its restriction to (0,t) belongs to llZ a (L°°(0,t)) for every t and 



A Q (ff): = — r SUp \(DlZ^)( S )\. 

1 (1 - a) 0<s <t<T 



1 , 

<:p— wT S l- a ,oo,T < oo 



We also denote by W^' (0,T) the space of measurable functions f on [0,T] such that: 



ds + / —, -—f-dyds < oo 



lo * a Jo Jo (s~y) a+1 
The restriction of / € Wq A (0,T) to (0,i) belongs to ig^L^O,*)) for all t. 

2.2 Taylor expansion of the solution 

For the convenience of the reader, we first recall the main result concerning the existence and 
uniqueness results for solutions of the differential equation driven by Holder path with Holder 
exponent f3 > 1/2. Several authors made contributions to this existence and uniqueness result. 
It has been first proved by M.Zahle in [13]. D. Nualart and A.Rascanu obtained the same result 
independently in [8]. Also, a different approach was introduced by A. A. Ruzimaikina to prove this 
existence and uniqueness theorem in [10]. The following result is due to Nualart and Rascanu ([8]) 
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Theorem 2.1 Let < a < 1/2 be fixed. Let g G W 1 ~ a ' 00 (Q, T; R d ). Consider the deterministic 
differential equation on M. n : 



X+ — X, 



+ f b l {s,x s )ds + y] f a^(s,x s )dgi, t e [0,T] (2.2) 

JO ~T JO 

— > M are measurable functions 



3=1 



i = 1, • • • , n where xq G R n , and the coefficients <r M , 6* : [0, T] x 
satisfying the following assumptions with p = l/a,0<r],5<l and 



■ 1 s 

< a < a = mm <j -,r/, — — 



1. o~(t,x) = (cr 1 ' 1 (t,x)) nX d is differentiate in x, and there exist some constants < n, 5 < 1 and 
for every N > there exists Mjv > suc/i i/iai i/ie following properties hold: 

'||<r(*,x) -a(t,y)\\ < M \\x -y\\, x G M n ,Vt G [0,T] 



\d Xi a{t,x)-d yi a{t,y)\\ < M N \\x - y\\ 5 \\x\\,\\y\\ < N,Vt G [0,T] 



Jk(t,x) -<7(t,i/)|| + \\d Xi a(t,x) -d yi a(t,y)\\ < Mo\\t — s\\ v ,\/t, s G [0, T] 



(2.3) 



There exists bo G Z/(0, T;R n ), where p > 2, and for every N > i/iere exists Ln > suc/i 
i/iai £/ie following properties hold: 



'\\b(t,x)-b(t,y)\\ <L N \\x-y\\, V\\x\\,\\y\\ <N,\fte [0, T] 
x)\\ < L \\x\\ +b (t), Vx G M. d ,Vt G [0,T] 



(2.4) 



Tfoen the differential equation \2. fy) has a unique solution x G W o '°°(0,T;W'). Moreover, the 
solution x is (1 — a) -Holder continuous. 

With this existence and uniqueness result in hands, we now come back to the differential equation 
we are concerned with: 



< 



dX t = Y,Vi(Xt)dyl 

8=0 

{X = x 



(2.5) 



Throughout the section, we make the following assumptions: 



1. The vector fields V-s are Lipschitz, C°°, with Lipschitz derivatives, and analytic on the set 
{x : \\x — xq\\ < C} for some C > 
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2. The driving path y : R> ->■ R rf+1 is /3-H61der continuous with (3 > 1/2. 

Clearly, from by Theorem 1, under our assumptions, this differential equation admits a unique 
solution (X t )t>o- 

Let / : R n — > R be a C°° function. By the change of variable formula, we have 



f(X t ) = f(x ) + E / VJ(X s )dyl 

JO 



Now, a new application of the change of variable formula to Vif(X s ) leads to 

f(X t ) = f(xo) + Y,Vif(x ) dyl+Y, / VjVifiX^dyidyi 
i=o Jo i,j=0 Jo Jo 

We can continue this procedure to get after N steps 

N 

f(X t ) = f( XQ ) + E E (^u ^/)(*o) / V + R N (t) 

for some remainder term Rjy(t), where we used the notations: 

1. A k [0, t] = {{h, ■■■,**)€ [0, t] k , < ti < t 2 ■ ■ ■ < t k < t} 

2. If / = (ii, ■ ■ ■ , ifc) G {0, 1, • • • , <i} fc is a word with length k, 

[ dy J =\ dy\\...dy\l. 

JA*[0,t] J0<t 1 <t 2 <--<t k <t 

If we let N — > +oo, we are led to the formal expansion formula: 

oo „ 

f(x t ) = f( Xo ) + E E Mi ■ ■ ■ v iM x °) / dy 1 

Now let us denote 7r J (x) = Xj, the j-th projection map. By using the previous expansion with 
/ = it 3 we get 

oo „ 

X l = 4 + E E • • • ^)(*o) / dy 1 



+ 0O 



= (U,— ,»fc) 

Therefore we have formally, 

X t = x + Ef fe W 
fc=i 

where 



u\=k jAk ^ 

This leads to the following definition: 



Definition 2.2 The Taylor expansion associated with the differential equation \2. 5|) is defined as 

oo 



k=l 

where 



Of course, at that point, the Taylor expansion is only a formal object in the sense that the conver- 
gence questions are not addressed yet. Our goal will be to provide natural assumptions ensuring the 
convergence of the expansion. Let us observe that the expansion can be written in a more compact 
way by using labeled rooted trees. For details on the combinatorics associated to these trees and 
the compact form of the expansion, we refer the interested reader to [7J. 

Another form of the expansion, which will be convenient for us, relies on the Taylor expansion of 
the vector fields and can be obtained as described below. We first need to introduce the following 
notations. Let bf (x — xo) a be the formal Taylor series of V\ at xq with the usual convention 

a£N n 

that x a = x" 1 • • • x" n if a = (ax, • • • , a n ) and x = (xi, • • • , x n ). Let now Cf 1 (hi, • • • , h k ) be the 
coefficient of y m in ^ bf (hiy H h h k y k Y . 

a£N n 

Consider the following rescaled differential equation, which depends on the parameter e: 



dXt = J2eV i (X$)dyt ^ 

oo 

Assume now that X\ admits an expansion in powers of e as xq + hk(t)e k for some h^s. Heuristi- 

k=l 

/ oo \ 

cally, because of the analyticity of the V(s, we can expand Vi I Xq + hk(t)e k at xq. Therefore, 
on the right hand side of (|2.6|) . we have: 



k=l 



eVi(Xl)dyi = * (Vi(x ) + C}(h{)e + Cf (h u h 2 ) e 2 + ■ ■ • ) dy\ 

8=0 1=0 



^2v i (x )dyt\e+ fecKhiWt) ^ + ' ' ' ) 

\i=0 J \i=0 / / 



Therefore, by identifying this expression with the left side of the equation we are let with the 
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conclusion that the hk(t)'s need to satisfy the following inductive system of equations: 

( d 
dh x {t) = J2Vi{xo)dyt 

d 

dh 2 (t) = Y J CHhi(t))dyl 



i=0 



(2.7) 



dh k (t) = Ct 1 (hi®, ■■■ , dy\ 



i=0 

oo 

If we let e = 1, we therefore see that xq + h k (t) formally solves the equation (|2.5p provided that 

fe=l 

the h k s satisfy the above inductive system. 
We have the following result: 

Proposition 2.3 Define h k inductively by the system \2. 7\ ), then for every k 6 N and i > me 
/iaue /i fc (t) = g k (t). 

Proof. 

First, let us observe that for every t > 0, e — ¥ X\ is a smooth function, indeed the equation 



i=0 



^0 = x o- 



can rewritten as 



dY t = Y J Vi(Y t )dy 



i=0 



Y = (e,x ), 

here Yt = (e, Xf) is now valued in Ex M n and V^(e, x) = (0, tVi{x)). The problem of the smoothness 
of X^ with respect to e is thus reduced to the problem of the smoothness of the solution of an 
equation driven by Holder paths with respect to its initial condition. Theorems that apply in this 
situation may be found in [6] and [9]. 

From the definition of hk(t) which is given by the system (j2.7|) . it is then clear by induction that 



h k (t) 



1 d k Xl 



(0). 



k\ de k 

On the other hand, as we have seen before, the repeated application of the change of variable 
formula to the smooth function f(x) = 7r J (x), gives 

N 



X 



£J _ rrJ 



l + E^ E (Vh-Vi^Xxo) [ dyt + 



JV+1 oe,j 
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where 



J=(ll,— ,Jjv+l) 1 J 

Therefore we have 

JV 

X{ = x + Y,9k(t)e k + e N+1 R%(t)- 

k=l 

d^ X € 

Computing , k * (0) by using the previous expression immediately gives /ifc(t) = □ 
2.3 Convergence of the Taylor expansion 

We now address the convergence questions related to the Taylor expansion. We use the same 
assumptions and the same notations as in the previous section. In particular we remind that the 
vector fields V-s are Lipschitz, C°° with Lipschitz derivatives, and analytic on the set {x : \\x—xq\\ < 
C} for some C > 0. 

2.3.1 A general result 

We first have the following general result. 

Theorem 2.4 Let xq + Y^k=l9k{p) be the Taylor expansion associated with the equation \2. 5|) . 
There exists T > such that for <t < T the series 

oo 

Z>*(*)ii 

k=l 

is convergent and 

oo 

X t = x + £flfc(*)- 

k=l 

Proof. 

Let us fix p > 0. For zdC, \z\ < p, we consider the complex differential equation 



dXt = Y J zV l {Xt)dy\ 

i=0 

which, by analycity of the V^'s, is well defined at least up to the positive time 

T{p)= inf inf{t > 0,Xf <£ B(x ,C)}. 

z,\z\<p 

As we have seen before, for a fixed t > 0, the map z — > X| is seen to be smooth. We claim it is 
even analytic. Indeed, differentiating with respect to z the integral expression 

,t d 



Xt=x + [ J2zVi(X z s )dyi, 
Jo 7T n 



8=0 

9 



immediately gives -gJ- = 0. That is z — > X* is analytic on the disc \z\ < p. 

On the other hand, xq + Y1T=i zk 9k(t) is precisely the Taylor expansion of Xf. Choosing then p > 1, 
finishes the proof. □ 

2.3.2 Quantitative bounds 

Theorem 12.41 is very general but gives few informations concerning the radius of convergence or the 
speed of convergence of the Taylor expansion 

oo 
k=l 

In this section, we shall be interested in more quantitative bounds. 

oo 

Lemma 2.5 Let x$ + <?/■(£) be the Taylor expansion associated with the equation \2. 5|) . For 

fc=i 

r > 1, we define 



T c (r) = M \t, J2 rk W9k(t)\\ >Cj . 

For every t < ?b(r) we have 

oo 

X t = x + ^2g k (t). 
fc=l 

Remark 2.6 By using a similar argument as in the proof of Theorem \2.4l it is seen that in the 
previous theorem, we always have Tc(r) > 0. 

Proof. Consider the following functions: 

Gi(e) = Vi ^ + f> fc ^(£)^ 

Then by the definition of Tc(r), if £ < 7b(r), Gi{e) is analytic on the disc {e € C : ||e|| < r}. 
Hence for any r' < r and if ||e|| < r', we have: 

171 — 1 

Gi(e) - (V-(xq) + ^( 3 i(£),52(£) ; <7 fe (£)))e A 



k=l 

oo 



(2i 



k=m 

<K\R m (2.10) 

where M % r , = sup Gi(e) < sup \\Vi(x + xq)\\ = K % And R m denotes the tail of the convergent 
||e]|<r' IN<C 
oo 



geometric series (~) ' ^ e know that — * 0, as m — * +oo. Now by the definition of gk{t) 

k=l 

we write: 
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d m 



k=l i=0 k=l J ° 

d m „t 

= E e E / (^ 1 Ct 1 (9i( S ),92( S ),.:,g k -i(s)))dyl 

i=0 k=l J ° 

d „ t I m-1 \ 

= E e / [Vi(x Q )+^SC}fa(8) > ....,g k (s)))dyi 

i=0 J ° \ k=l J 

oo 

Since when t < Tc(r), the series e k gk(t) < C > w e replace t by t f\Tc{r) in the previous equation: 



fe=i 



m rf ftAT c (r) ( m- 1 \ 

Y,e k 9k(tAT c (r)) = J2t / ^(xq) + ^ e fc Cf .... )5fc ( S )) ] ^ (2.11) 

fc=l 8=0 ^° V 



fc=l 



with the estimate (|2.10p . we have when e < r', 



TO— 1 



Gi (e) - K-(x ) + E C ^ (9iV),92{t), ....g k {t)) e k 



k=l 



< K l R ri 



(2.12) 



Therefore, since r > 1 we can choose 1 < r' < r and let m — > oo, then the right hand side of (|2.12p 
tends to 0. Therefore, by letting m — > oo in equation ([2. lip we have: 

00 d rtAT c (r) 

x + t k 9k(t A T c {r)) =x + Y *Gi(t)dyi 
fc=l i=0 ^° 

d ftAT c (r) / oo \ 

=2; o + E / eFi x ° + E d f* 

i=0 ^° V fc=l / 



And the result follows when we take e = 1 



□ 



The crucial estimates we need to control the speed of convergence of the Taylor expansion are 
provided by the following proposition. 

Proposition 2.7 Let a such that 1 — (3 < a < h and fix T > 0. For < t < T, we have 



L 



dyll • • • dyl n 



0<ti<...<t n <t 



< 



T(l-2a) ,- r ,n-l. , .... .. .,.„-.! 



r(n(l - 2a)) 



C a (T)"- i A a (r,y) ft - 1 [|y|| a>Tl< 



where 



Aa(T,y) := 



C«(T) 



1 



r(l - a) < s <t<T 



sup ||(A- Q yt-)( S 



i 



+ T a (1 +T a )T(l - 2a)T 



l-2a 



a(l — a) 

0<t<T V JO {t~ S) + 
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Proof. Since Holder path with exponent (3 > 1/2, and since a + (3 > 1, we have 

y\ 6 Ty^" Q '°°(0,r) for * = 0,1,2- •• ,d. For t < T, we denote I fc (t) = / dy\\- ■ ■ dy\ k k , 

Jo<t 1 <...<t k <t 

k = 1, • • ■ , n. Clearly, 4(i) G W^'^O, T) for = 1, 2, • • • , n. Therefore, Proposition 4.1 in Nualart- 
Rascanu ([8]) gives 

i4 + iwi+y o (f _ a)1+a & 

<A a (T, y) + T«) jf ((t - + rD (|4(„)| + f ^T#^) ** 

Using now the basic estimate 

(t - r)~ 2a + r~ a < (1 + T a ) t 2a (t - r)- 2a r- 2a , 

we get 

Iterating this inequality A; times, we get: 
' |4+ i (*) - h +i(s) 

(t- s y 

< Aa(T, y) k (^^y + 7°) (1 + T a ) k \\y\\ a , Tj00 J* (t 2a (t - r k )- 2a rf a ) jf' (r 2a (r k - n^y 2 ^ 2 *) 

[ k ~ l (rl a -i(r k -i ~ r k - 2 )- 2a r~ 2a 2 ) ■■■ T {r 2a (r 2 - n)' 2 ^) d n dr 2 ■ ■ ■ dr k 
Jo Jo 

(i \ k i't rrk r r k—i 
+ t«) (1 + T a ) k \\y\\ a , T ,oo / {t 2a (t - r k )~ 2a ) / (r k - r k ^y 2a / (r fc _j - r fe _ 2 )- 
"(1 -a) J Jo Jo Jo 



2a 



'■■2 



(t2 — ri) 2a r x 2a dr\dr 2 ■ ■ ■ dr k 



o 



ri 

Note that by the substitution u = — , we have: 

r 2 



f 2 (r 2 - n)- 2a r^ 2a dn = r\- 4a f (1 - u)" 2a u~ 2a du = r£ -4a B(l - 2a, 1 - 2a) 

JO JO 



where 

filer. ii\ = 

T{x + y) 



B(x,y) = ^4= I (l-^-V-^u 



o 



12 



is the Beta function. 

Then plug r\~ Aa B(l — 2a, 1 — 2a) into the previous integrand and iterate this process, we get: 
\I k+1 (t) - I k+1 (s)\ ^ 

l4 + l(i)l+/ Q {F ^y T - a ds 

< A a (T, y) k ( — - + T a ) (1 + T a ) k ||y lU^ooir 1 - 2 ")^ TT B(l - 2a, i(l - 2a)) 
Va(l-a) / a a 

^ AJT, yf (-* + T°Y (1 + ItolkT^tT-'-)' fl r( i:. 2 °!^" i 

Va(l-a) / ax r((t + l)(l-2a)) 



' ' ' 2 ° ' A a (T, y) fc f 1 + T a ) " (1 + T a ) k T(l - 2a) fe (T 1 - 2Q ) fc ||y|| Q , T , c 

\a(l — a) J 



r(fc + 1)(1 - 2a)' 

□ 

We are now in position to give an estimate of the speed of convergence of the stochastic Taylor 
expansion. The following lemma is taylor made for that purpose: 

Lemma 2.8 Let a such that 1 — j3 < a < | and 7 < 1 — 2a. There exists a constant K a>1 > such 
that for every N > and x > 0, 



r(fcr) 

; a)) 

^ r((l-2a- 7 )iV) ' 



^ T(k(l-2a)f ~ Ka ^\ 

k=N+l y v " 



e 2xT= ^, ifN = 



jsj o-r 1 — ' 2a — 7 

* e i/iV > 1 



Proof. We make the proof for iV > 1 and let the reader adapt the argument when N = 0. We have 

x N ^ r((jfc + N + l)7)r((l - 2a - 7 )iV) ^ 



r( 



T((l - 2a - 7 )JV) ^ r((fc + iV + 1)(1 - 2a)) 

- Ca '"T((l-2a- 1 )N) £ r((fc + l)(l-2a- 7 )) 
From Lemma 7.7 in [8], we have for every x > 0. 

00 ^fc /| c 2 1 

r((l + fc)(l-2a-7)) ~ 1 - 2a -7 6 

This concludes the proof. □ 

We can now give convenient assumptions to decide the speed of convergence of the Taylor expansion. 
We remind the reader that if I = (ii, ■ ■ ■ , i k ), we denote 

Pi = (V tl ---V lk ^)(x ). 

We also use the notations introduced in Proposition 12.71 
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Theorem 2.9 Let a such that 1 — j3 < a < \. Let us assume that there exist M > and 
< 7 < 1 — 2a such that for every word I E {0, • • • , d} k 

\\Pi\\ < r(7|/|)M |71 . (2.13) 

Then, for every r > 1, Tc(r) > and there exists a constant K a ^ t M,d > depending only on the 
subscript variables such that for every < t < Tc(r) and N > 1, 



/ N N 

V k=l / 



< K tf I, ii (M(d + l)A Q (t, y)C a (t)) N 2 (M(d+i)A a {t,y)C a (t))^ 
S n a , liM , d \\y\\a,t,oo r((l - 2a - 7 )JV) 



In particular for every < t < Tc(r) 



+oo 



-Xt = x + E#fc(*)- 



(2.14) 



(2.15) 



k=l 



Moreover if C = +oo, that is if the Vi 's are analytic on 1", then \2.1J^ and 12.15]) hold for every 
t > 0. 

Proof. Let us fix T > and r > 1. For every t G [0, T], by Proposition 12.71 we have: 

oo 

E rfc i^(*)ii 

dy 1 



k=l 



k=l \I\=k 
oo 



A fe [0,t] 



< J2 r k (d + l) k M k T(k 7 )A a (T, y)^ 1 ^^)^ 1 || y || a>T>0O ^ 



r(l - 2a) 



k=l 



(fc(l -2a)) 



< ||y||a,r,oor(d + 1)MT(1 - 2a) £ (r(d + l)C a (T)A a (T, y)) 



fe-i 



r(fc 7 ) 



fe=i 

oo 



r(Jfe(l-2a)) 



|y||a,T,oo^(d + l)AfT(l - 2a) 2^ ^ttt; — Tv; B(kj, k(l -2a- 7)) 



k=l 



r(Jfe(l - 2a - 7)) 



< 5(7, 1 - 2a - 7)||y||a,T,oor(d + 1)MT(1 - 2a) g 



fe=i 



(r(d+l)C a (r)A a (T,y) ) 
r(fc(l-2a-7)) 



fe-i 



Now note that by our assumption, 1 — 2a — 7 > 0, then the series on the right hand side of our 
inequality converges. This implies that 



T c (r)=M[t, J2 rk \\9k(t)\\>C) >0. 



k=l 



We conclude from Lemma 12.91 that for every < t < Tc(r), 

+00 

x t = x + e 



fc=i 
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If moreover C = +00, that is the radius of convergence of the Taylor expansion of the V^s is +00, 
the previous equality holds for every t > because, limc , _ > .+ 00 Tc(r) = +00. 
We then get for < t < Tc(r), thanks to Lemma 12.81 



k=N+l 



£ E E i*\ 

k=N+l \I\=k 



A k [0,t] 



dy 1 



< f; (d + l) fe M fc r(A: 7 )A Q (t,y)^ 1 C Q (T) fc - 1 ||y|U, 4 , oo! i^-^ )) 



k=N+l 



<K ii-.m (M(d + 1) A Q (t, j/)C a (t))^ 2 (Af (d+l)A n fa/)Cq («)) ^ 



□ 



2.4 The Lie group example 

In this section, we will see an application of our main results in the context of Lie group. Let us 
denote by G a connected Lie group. For any g € G, we have left and right translation L g and R g , 
which are diffeomorphisms on G, defined by L g (h) = gh and R g (h) = hg for any h £ G. We say 
that a vector field V is left-invariant on G if for any g £ G, 



(L g UV) = V 

Let us denote by q the Lie algebra of G, that is the set of left-invariant vector fields on G and 
by T e G the tangent space of G at identity. We know that the map V —> Vie) provides a linear 
isomorphism between g and T e G. Now let Vi, ■ ■ ■ , Vd € JJ be left-invariant vector fields and let us 
consider the following differential equation on Lie group G, 

,«. = E*<*.» (2 . 16) 

X = e 

where y : M>o — > M. d+1 is a /3-H61der path with ft > 1/2. The vector fields are not supposed to be 
globally Lipschitz, so that we are not exactly in the previous framework. However, the vector fields 
Vi's are smooth and analytic and the previous equation always admits a solution up to a possible 
positive explosion time r. 
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Remark 2.10 If&= GL(n,M), we know the associated Lie algebra q = M(n,M). In this case, the 
Taylor expansion associated with \2. 1 6\) is convergent for every t > 0. We actually have 7 = in 
the convergence criterion of Theorem \2.9[ 



Before we state our main theorem in that context, let us introduce some notations (see [2] and [1] 
for further details). 

• If / = (ii, • • • , ik) we denote by Vj the Lie commutator defined by 

V I = [V h ,[V i2 ,...,[V ik _ 1 ,V ik }...}. 

• 5^ is the set of the permutations of {0, k}; 

• If a € Sk, e(ci) is the cardinality of the set 

{j e{l,...,k-l},a(j) >a(j + l)}, 



°es k k 2 ) J o<ti<...<t k <t 



e{a) 

The main theorem of this section is the following: 
Theorem 2.11 There exists T > such that forO<t<T, 

(00 
£ £ Aj(y) t V^ ) , 
fc=l (n,...ifc) 



is 



the solution of the equation \2.1b\ 



Proof. Let us denote by ad the map (adx) (z) = [x, z\. Consider the set g e = {x G T e G , s(x) fl 

2itK C {0}}, where s{x) denotes the spectrum of adx. It is known that g e is an open set 

( I — e~ adx \ ( I — e~ adx \ 

of T e G and Vi £ j e , I I is invertible, where I I is defined as power series 

\ adx J \ adx / 

y — — (— adx) . Now we define vector fields U% on g e by UAx) = ( , I (VAe)) where 

£^ (k + iy. \I-e adx j 

adx \ , , . „ (I — e ~ adx 



denotes the inverse of 



I — e~ adx J V adx 

We claim that for Vx E g e we have: 

(T x exp)[Ui{x)} = Vi(expx) 
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In fact, since V(s are left-invariant vector field on G, we have 



(I - e ~ adx \ 

(T x exp) [Ui(x)\ =(L(expx))* o I — —— — j [Ui(x)\ 

= (L(expx)),(F,(e)) 
=\^(expx) 

Now if we consider the differential equation on g e : 



dY t = Y j U i {Y t )dy\ 



k y = o 



(2.17) 



we know that Ui(x) are analytic on Q e . Therefore from Theorem 12.44 there exists a positive time T 
such that for < t < T, 

oo 

Y t = Y J 9k{t), 

k=l 

where YlkLi 9k{t) is the Taylor expansion of the equation defining Y t . It is proved as in [2] and [TT] 
that 

g k (t)= Mv)tVi. 
(n,...ifc) 

Now we claim that by the chain rule, we have Xt = expiYt) when t < T. In fact: 

dX t =(T Yt exp)dY t 
d 

=(T Yt exp)[J2Ui(Y t ))dyi 

i=i 



= Y J T Yt exp[U l (Y t )}dyl 

i=0 
d 

= X)^(exp(y t ))dyj 
i=0 
d 

i=0 

with Xq = exp(Yo) = e 



□ 
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Remark 2.12 It seems difficult to obtain precise universal bounds on the radius of convergence of 
Sfc=i i k ) h-i (y)tV[ a s it is the case for the radius of convergence of the Magnus expansion of 

absolutely continuous paths. Indeed, rough estimates lead to 



e(a) 

for some constant C > 0. We get therefore the bound 



\Ai(y)t\ < ^ C^tf-'A^y)^ 1 ] 



v\ 



0,00 1 



2 k ' r(fc(l - 2a)) 

which unfortunately gives a convergence for the series for a = 0, which corresponds to the case 
where y is absolutely continuous. 

3 Taylor expansion for stochastic differential equations driven by 
fractional Brownian motions 

We now turn to the second part of this work and applies our result to stochastic differential equations 
driven by factional Brownian motion paths. 

Let us recall that a d-dimensional fractional Brownian motion with Hurst parameter HE (0, 1) is 
a Gaussian process 

B t = {Bl--- ,Bf), t>0, 

where B 1 , ■ ■ ■ ,B d are d independent centered Gaussian processes with covariance function 

R{t , s ) = ±{s 2H + i 2H -\t-s\ 2H ). 

It can be shown that such a process admits a continuous version whose paths have finite p variation 
for 1/p < H. Let us observe that for H = ^, B is a standard Brownian motion. 
In this part, we will consider the following stochastic differential equation: 

d 

dX t = y^Vi (X t ) dBi 

^ v ' 1 3.18 

Xq = x 

where the V^'s are Lipschitz C°° vector fields on R n with Lipschitz derivatives and Bt is a d- 
dimensional fractional Brownian motion with Hurst parameter H > 1/2 with the convention that 
B = t. 

Let us remember that when H > 1/2, Bt is almost surely Holder path with Holder exponent strictly 
less than H. In particular, we can pick 1/2 < (3 < H such that Bt is (3 Holder path. Hence, we can 
define the Stochastic Taylor expansion associated with (|3,18p just as in the deterministic case by 



k=l 
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where ^(t) = £ P\ [ dB 1 , P\ = (V H ■ ■ ■ V ik rf)(x ) and J Ak dB 1 = [ dB\\ ■ 

| 7 | =fc J A k [0,t] ' J0<ti<t 2 ---<t fc <t 

As in the previous section, we assume throughout this section assume that the V^s are Lipschitz 
C°° vector fields on M. n with Lipschitz derivatives and analytic on the set {x : \\x — xq\\ < C} for 
some C > 0. 

Now let us state the following convergence theorem for solution to the stochastic differential equation 
(I3.18P , which is an immediate consequence of Theorem 12.41 

Theorem 3.1 There exists an almost surely positive stopping time T, such that for < t < T the 

series 

oo 

Ikoll +^11^)11 

i=l 

almost surely converges and almost surely 

+ 0O 

X t = x + ^2g k (t). 

k=l 

3.1 Quantitative deterministic bounds 

In this part, we will derive a convergence criterion for the stochastic Taylor expansion associated 
with the SDE (|3. 18|) . The next theorem is just an direct application of the deterministic theorem 



Theorem 3.2 Let a such that 1 — H < a < \. Let us assume that there exist M > and 
< 7 < 1 — 2a such that for every word I G {0, • • • , d} k 

\\PjW < r( T |/|)ikri 7 i. (3.19) 

Then, for every r > 1, Tc(r) > almost surely, and there exists a constant K a ^ t M,d > depending 
only on the subscript variables such that for every < t < Tc{r) and N > 1, 

Xt- (so + |>(t)) <K a ,,, M AB\\ a ,t J M{ ^ 

(3.20) 

In particular for almost every < t < Tcif), we have almost surely 

+ 0O 

X t = x + Y J 9k(t). (3.21) 

k=l 

Moreover if C = +oo, that is if the Vi 's are analytic on 1", then 113. 20\) and 113.21]) hold for every 
t > 0. 
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We can observe that the error estimate provided by the previous theorem is in L p for every p > 1. 
Indeed, Lemma 7.4 in Nualart-Rascanu ([8]) says for any (3 < H and T > there exists a positive 
random variable rjp t T such that E(\r}p t T\ p ) < oo for all p > 1 and for all s,t G [0,T] 



|-B(t) -B(s)\\ < r]p iT \t-sf a.s. 



We now have, 



l-S|loToo= SUP I ll-Btll + / -7 — TTT— 

< sup (vp,T^ + [ rip^t-sf-^ds] 

0<t<T V JO / 

So, E||S||^ Too < oo follows from the fact that E(\ i qpp\ p ) < 00 f° r all p > 1 and for all s,t G [0,T]. 
As for A a (T,B), we have the following estimate 

A a (r,B)= T , /i 1 — sup HD^Bt-WII 



r(l - a) <s<t<T 

1 / 1 1 Bt — B s || . /" || S M — B s 

r(i - a)r(a) <!<?<T U-*) 1 -" + (1 " a) X ( U - S )2- 



< — sup (npr(t - s)^*- 1 + (1 - a)r,p >T [ (u - sf +a ' 2 du) 

Cs<t<T V J s J 



1 

T(l-a)r(a) . , 



Sup \(t-sf +a - l + 7T—{t-s) 

r(l - a)r(a) < s <t<T V a + /3-l V ; 



< 



/3T' 



Vp,T 



(a + /3 - l)r(l - a)r(a) 

Note that by assumption we know a + /? — 1 > and E|A a (i?)| p < oo follows by the fact that T)p T 
has a Gaussian tail from Fernique's theorem. 

Remark 3.3 Note that in this argument we didn't include the drift term when we showed the LP 
integrability. However, the argument can be modified trivially if we include drift term, since the 
fractional derivative oft with order 1 — a > is bounded in any closed interval [0,T] 

3.2 Quantitative probabilistic bounds, L 2 error estimate 

In this part, we will develop a convergence criterion for the stochastic Taylor expansion which is 
based on the L2 estimates of the iterated integrals of fractional Brownian motion. Such estimates 
may be found in Proposition 4.8 in [7j. They rely on the computation of E ^/^mp 1 ] dB 1 ^ established 
in [3] and the isometry property of the divergence operator. 
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Compared with the deterministic approach, the probabilistic method improves the speed of conver- 
gence in I? when H G (|, |). 



For convenience in this section we restrict ourselves to the case where there is no drift, that is we 
assume Vq = 0. By using heavier notations, the arguments we use may be extended when Vq 7^ 0. 

Let us first remind the result by Neuenkirch-Nourdin-Tindel in [7]. 

Proposition 3.4 Assume that H > 1/2 and let an = H(2H — 1) and K = yjj^- we have the 
following inequality: 









E I 


[ dB 1 






JA m [0,t] 





ml 

It is interesting to observe that this decay in m is faster than the decay obtained by the deterministic 

'1 3^ 
v2> 4J 



method when H £ §). We now prove: 



Theorem 3.5 Let us assume that there exists M > andO < 7 < 1/2 such that \\Pi\\ < M^'d/I!) 7 . 
Then for every r > 1, Tc(r) > almost surely and for < t < Tc(r) ; we have in L? and almost 
surely 

+00 

X t = x + ^2g k (t). 



k=l 



Moreover there exists a constant C 7 > such that, 



E 



N 



k=l 



where 



x* 



k=0 



(fc!)V2- 



Moreover if C = +00, that is if the Vi 's are analytic on M n , then we may take Tc{r) = +00 in the 
above inequality. 
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Proof. Let p > 0. We have: 
Pi 



E 

.\I\=m 



dB 1 



A m [0,i] 



E 

J/|=m 



Pi 



/ dB 1 

JA m [0,t] 



2-i 



< dm V E 

H \I\=m 



L 



P 7 / dB 1 

A m [0,t] 



\I\=m 



dB 1 



A m [0,t] 



d 



<(- ) E wpiW 2 ^^ 

\P 2 J ^ m rn\ 

<{-) M^j^—t^* 
\ p 2 / (m!) 1 ^ 



Therefore: E P E 

rn=l \\I\=m 



Pi I dB 1 

A m [0,t] 



> p m \ < oo. From Borel Cantelli Lemma, we deduce 



that for every r > 1, ^ r m ||g m (i)|| converges uniformly on any interval [0,T] almost surely. From 

m=l 

the results of the previous section, it thus implies that for every r > 1, Tc(r) > and that for every 
< t < Tc(r), we have almost surely 



Finally, we have 

X t - ix + ^2g k (t) 



E 



N 



k=l 



k=l 



1/2 



K k 



,t<T c (r)) < £ ^^pt^M'W 

k=N+l \I\=k * 



+oo 



- S (^1)1/2 

k=N+l v > 



(dKt 2H M) k 



< a 



{dMKt 



7 

2H\N+l 



%(dKMt 2H ). 



((N + 1)!) 1 / 2_ t 



□ 
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